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The cold(neutronless) fission of ^^^Cf is studied in the frame of a molecular model 
in which the scission configuration is described by two aligned fragments interacting 
C3 . by means of Coulomb ( + nuclear) forces. The study is carried out for different dis- 

^^ \ tances between the fragments tips and excitation energies. For a given deformation, 

^^ ' the fragment's total energy is computed via the constrained Hartree-Fock + BCS 

0> . formalism. The total excitation energy present in the fragments is supposed to con- 

tribute only to the fragments deformation and the asymptotic value of the kinetic 
energy is equated to the inter-fragment potential at scission. These two constraints 
yield not more than one or two fission channels for a fixed tip distance and excitation 
O , energy. Discarding those fission channels corresponding to a disequilibrated sharing 

of the excitation energy between the two fragments, we were able to establish the 
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most likely scission configurations for a specified excitation energy. 
PACS numbers 21.60 Gx, 21.60. Jz, 24.75.+i, 25.85. Ca 
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1 Introduction 



In last time a renewed interest in the spontaneous fission (sf) of ^^^Cf arised in 
connection with modern experimental techniques, based on large Ge detector 
arrays, which allow a better determination of the mass, charge and angular 
momentum content of the fragments [1]. 

A particular attention has been payed to the limiting case of cold fission, when 
no neutrons are emitted and the energy released in the reaction is converted 
entirely in the kinetic energy of the fragments. Some features of this process 
have been very recently explained satisfactory using cluster like models[2,3]. 
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In these models it is assumed that at scission the fragments have very compact 
shapes, close to the ground state and thus they are carrying very small excita- 
tion energy. The scission configuration consists of two co-axial fragments with 
a certain distance d between their tips. In the model proposed by Gonnenwein 
et al.[4], the cold fission was studied by determining the distance dmin of the 
closest approach between the two fragments, when the Q-value equals the in- 
teraction energy. This model predicted the smallest tip distance (bellow 3 fm) 
for fragments, with mass numbers between 138 and 158 and around the double 
magic ^^^Sn, emerging in the sf of ^^^Cf. Small tip distances were interpreted 
as a sign of cold fission due to the higher interaction energies at scission. 

In the past the scission-point model succeded also to explain roughly some 
basic observables of low-energy fission. Based on the assumption of statisti- 
cal equilibrium among the collective degrees of freedom at the scission point, 
Wilkins et al. [5] calculated the relative probabilities of formation of comple- 
mentary fission fragment pairs from the relative potential energies of a system 
of two coaxial, quadrupole deformed liquid drops, with shell corrections taken 
into account. The distance between their tips, the intrinsic excitation energy 
and collective temperature were choosen as the free parameters of the model. 
In this way the general features of the distributions of mass, nuclear charge 
and kinetic energy in the fission of various nuclides, ranging from Po to Fm 
were well reproduced. Using similar ideas, Norenberg[6] computed the level 
schemes, equilibrium deformations of the fragments, total energies and charge 
distributions of ^^^U, 240-242p^ using the BCS wave-function in the description 
of the ground state. 

In this paper we generalize the static scission-point concept of nuclear fission 
model in such a way that instead of describing the fragments as two deformed 
nearly touching liquid drops with shell corrections taken into account, we in- 
corporate the fragments shell structure by means of the self-consistent Hartree- 
Fock(HF) method. For the given binary splitting ^^^Cf^^^'^Mo + ^"^^Ba we first 
established the equilibrium deformations of the two fragments by seeking the 
HF minimum and next their total energy for various deformations is computed 
by constraining their quadrupole moments. The two fragments are considered 
as coaxial with distance d between their tips. 

One of the basic approximation employed in this paper was that the inter- 
action energy at scission is transformed into kinetic energy of the fragments 
at infinity. Thus, all the excitation energy present in the fissioning system is 
accounted by the deformation energy. This amounts to neglect that part of 
the energy released at the descent from saddle to scission which is spent on 
heat. Thus, our study concerns mainly the low-energy domain of sf including 
the limiting case of cold fission [4,7]. 

By using the above mentioned constraints we were able to deduce the possible 



shapes of the fragments for various tip distances and total excitation energies 

E*. 



2 Molecular Model of Low Energy Fission 

2.1 Energy Balance 



In the sf of ^^^Cf the fragments are born with a certain deformation and will 
carry a total excitation energy E*, gained during the descent from saddle to 
scission which will be dissipated by means of neutron and gamma emission [8] 

Q = Vsc^ss + TKEpre + E ^defi^ + E ^™*(^) (1) 

1,2 1,2 

where Vgdss = Vcoui + K.«d represents the fragments interaction energy at 
scission. For Vcoui we choose the form corresponding to two deformed homoge- 
nously charged nuclei with collinear symmetry axes with a distance R between 
their centers [9] 
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The variables xi and X2 are expressed in terms of the semiaxes a and b 

"1,2 "1,2 /o\ 

a;i,2 = -^^ (3) 

The above double series is converging for |a;i| + |a;2| < 1 and the final result is 
given, according to [10] : 
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ln(l + xi + 0:2) (1 — 3^1 — 3^2 + 12xia;2 — 4x^ — Axl) \ (4) 

For the attractive nuclear potential we choose the proximity formula for two 
nuclei with a finite surface thickness[ll] 

Vnuci = 47ri?7$(C) (5) 



The explanations of the different quantities entering in the above formula can 
be found in [12]. The prescission kinetic energy TKEpre is taken to be zero, 



an assumption which proved to be reasonable for low-energy fission[5]. Also, 
that part of the excitation energy which is transformed into internal excitation 
energy E*^^, is neglected. According to calculations based on time dependent 
quantum many-body calculations, the intrinsic excitation energy accounts for 
less than 15% of the collective energy gain in going from the saddle to the 
scission [13]. Schiitte and Wilets [14] gave also an upper bound for -E*„t) which 
is still small compared to the total excitation energy E*. 



2.2 The deformation energy in the frame of the Hartree-Fock + BCS method 



That part of the excitation energy which goes into the deformations of the 
fragments was denoted in eq.(l) by Ed^j- In the study of sf properties, E^ef 
is taken usually as a sum of the liquid drop model (LDM) energy, and the 
shell and pairing corrections [15]. In this paper the deformation energy Edef 
of the fissioning system at scission is referred to the HF+BCS energy of the 
two fragments in their ground states 

Edef = Ehf+bcs{Ni, Zi, Pi) — Ehf+bcs{Ni, Zi, Pi ') 

+ Ehf+bcs{N2, Z2, P2) - Ehf+bcs{N2, Z2, PI-'-) (6) 

Obviously this is a more general approach. The LDM, which is based on a 
semiclassical description of the nuclei, supplemented by the shell-effect cor- 
rective energy, is only a poore substitute for a self-consistent calculation. One 
of the main advantages of the self-consistent HF+BCS calculation is that it 
provides simultanously both the single-particle and semiclassical properties of 
nuclei. The general properties of the Hartree-Fock method were reviewed in 
[16,17]. 

In our study for the HF part of the interaction we choosed the Skyrme inter- 
action SHI [18], which succeded to reproduce satisfactory the single-particle 
spectra. The difference between the binding energy computed with SIII and 
the experimental one appears to be, for a large number of nuclei, ~ 5 MeV 
[19]. It also produces a fairly well N — Z dependence of the binding energy [20]. 
The present study envisages nuclei that are not in a closed shell configuration. 
Thus, the level occupations will have a large effect on the solution of the 
HF equations. Usually the HF method is extended to the Hartree-Fock Bo- 
golyubov (HFB) formalism by using a mixture of different configurations in 
place of a single Slater determinant. However, when dealing with a Skyrme 
force which has been simplified such that the bulk properties of the nucleus 
are reproduced, one would have to introduce additional parameters in order 
to guarantee that sensible pairing matrix elements are obtained. 



Following Vautherin [21] we assign to each orbital (pk an occupation n^ 



•'fc; 



where ul + 



1, ui = Uk and v^ = —Vk- In terms of the density p(r) 



2Z]A:^fc|0A:(^)P the HF+BCS total energy, that has to be minimized reads 
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where 
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is the expectation value of the kinetic energy, V = tr{pv) enters as Hartree- 
Fock-like potential, v being the antisymmetrized effective two-body interac- 
tion. The primed sum J2' denotes a sum over all HF orbitals having projections 
of the total angular momentum j on the z-axis ^2^ > 0. To the total energy 
we added the pairing energy 
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(9) 



For BCS-like calculations, the matrix elements of v between HF states is taken 
to be constant 
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Varying the normalized single-particle wave functions (j)k and their amplitudes 
Vk under the additional constraint XrUk {^Tk,Tnk — Nr) , (r = p,n), ensuring 
that on the average the system contains the correct number of neutrons A^ 
and protons Z, we are lead to the standard HF and BCS equations [21]. 

The occupations Uk are determined at each step of the HF iterative calcula- 
tion using the HF eigenvalues Sk and they are employed at the next step to 
construct the HF field. The pairing force constant is 
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MeV (r = p, n) 
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The constant Gqt was adjusted in such a way to obtain the experimental 
pairing gap 



A^ = GJ2 UkVk 



:i2) 



In the deformed HF calculations one have to optimize the basis which is 
choosen to correspond to an axial symmetric deformed harmonic-oscillator 
with frequencies uj_[_ and Uz- Such a basis is characterized by the deforma- 
tion parameter q = uj±/uJz and harmonic oscillator length b = Jmcoo/h, with 
Uq = ujj_Uz. The basis is cut off after Nmax major shells, where Nmax=^0 or 
12 for the nuclei emerging in the sf of ^ssqj [22]. 

The next step consists in mapping out the potential energy curves by con- 
straining our HF+BCS calculations in which a quadratic constraint ^{Q—QoY 
is added to the energy functional (7) [23]. Here Qq is a specified value of the 
mass quadrupole moment. In Fig.l we represent the deformation energy curves 
of the nuclei ^°^Mo and ^^^Ba produced in the sf of ^^^Cf. 



3 Distribution of the excitation energy in the fission fragments 



The scope of this section is to seek the configuration of the system at scission 
for a fixed excitation energy E*. According to eq.(l), the interaction energy 
of two fragments with deformations Pi and P2 at scission is related to the 
excitation energy through the relation 

V{f3i,f3,,d) = Q-E* (13) 



where d is the tip distance and enters in the theory as a parameter. We equate 
this last quantity with the asymptotic kinetic energy TKE{oo). This relation 
is a conequence of the approximations that we made earlier, i.e. we neglected 
the prescission kinetic energy TKEpre and we forced all the available excitation 
energy to be stored into deformation 

E*(/3i,/32) = Edef iPi) + Edefm (14) 

where E^ef is computed according to eq.(7). Thus, for a given excitation energy 
we obtain two non-linear equations, i.e. eqs. (13) and (14). 

In fig. 2 we represented the excitation energy landscape (14), for the pair 
(^°^Mo, ^^^Ba). The deepest minimum corresponds to the prolate-prolate con- 
figuration (/^i, /32). At this point _E* = and fission proceeds by means of only 
one channel, customarly known as cold fission. This configuration has defor- 
mations /5("'^°^Mo)=0.370 and /3("'^'*^Ba)=0.270 which are very close to those 
computed in the frame of the finite-range droplet macroscopical model [24]. 
The non-linear equations, quoted above, admit this solution only for the tip 
distance d = 2.95 fm, a value very close to the border of 3 fm, allegated by 
the Tiibingen group, bellow which cold fission occurs [3]. 



When we increase the excitation energy, an infinity of solutions arise according 
to eq.(14). They have to be identified with the geometrical locus of points 
with equal excitation energy. However, the second constraint (13) is limiting 
drastically the number of {f3i,f32) pairs. In Fig. 3 we give the contour plots of 
the excitation energy and superposed on them the curves relating j3i to /32 for 
different tip distances. The intersection of such curves with the contour lines of 
equal excitation energy will give the physical solutions to our fission problem, 
i.e. for certain tip distance intervals, one get different scission configurations 
or channels. 

As one observe one get generally two solutions which are located mainly for 
low-excitation energy in the quadrant with /3i,/32 > 0. While for pure cold 
fission {E* = 0) one get a solution only for one d, when E* > one get solu- 
tions for several values of d. Naturally, one may ask next if all these solutions 
are likely to occur. For that one should look at the ratio of excitation energies 
between the two fragments. Calculations based on the cascade evaporation 
model predicted a ratio of the mean excitation energies E2/EI ^ 0.5 around 
the splitting 104/148 when approaching the limiting case of cold fission [25]. 
According to the same reference, disproportions in sharing of the excitation 
energy should be expected only in the vicinity of magic numbers. For our 
study case we are left only with few possibilities for a given excitation energy, 
which are listed in Table I for E*= 0, 2 ,4 and 6 MeV. 

In fig. 4 we give the fragments density contour lines for a fixed excitation 
energy, namely E*=2 MeV and different tip distances. 



4 Conclusions 



Based on a molecular model in which the scission configuration has to ful- 
fill two main energetic constraints, namely that the interaction between the 
fragments is converted totally into asymptotic kinetic energy and that the ex- 
citation energy of the fissioning system is accounted only by the deformation 
energy, we carried constrained HF+BCS calculations at zero temperature for 
the nuclei emerging in the low-energy fission reaction. For a fixed excitation 
energy we varied the distance between the tips of the fragments. Each case ad- 
mits not more than two solutions, i.e. to pairs of fragments deformations. The 
criteria which allowed us to select the valid scission configuration was that of 
the excitation energy distribution between the fragments. We discarded those 
configurations with a disproportionate ratio between the excitation energis 
of the two fragments. A carefull analyse exquibite roughly two regions of tip 
distance which can be assigned as valid scission configurations. The first one 
has the starting point aX d = 2.95 fm and goes up to 3.8 fm whereas the 
second is much narrow and is centered around d = 5fm. In this last case that 



should be less probable to occur, one of the fragments is emitted with oblate 
deformations. 

For excitation energies higher than those considered in this paper, the prop- 
erties of the nuclear system are described by a thermal average, the influence 
of the shell effects becoming thus less important. 

A limitation of the present approach is caused by the absence of higher mul- 
tipole deformations (octupole, hexadecupole, etc.). As have been shown very 
recently [2], the account of hexadecupole deformation provided the explana- 
tion of a whole region of cold fission for ^^^Cf. 

Also we intend to study the fragments angular momentum formation in these 
fragments based on a very recent proposal of some of us [26] together with the 
evolution of several collective variables during the post-scission motion. 
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Figure legends 



Fig.l The deformation energy curves of the nuclei ^^^Ba and ^°^Mo computed 
in the frame of the HF+BCS method with quadratic constrain for the mass 
quadrupole 

Fig. 2 Three-dimensional plot of the excitation energy E* for the pair (^^^Ba, 
^°^Mo) computed in the frame of the HF+BCS method. 

Fig.SGraphical solution of the non-linear equations (13) and (14). The inter- 
section of the solid curve with the contour lines provides two solutions in the 
particular case of the pair (^^^Ba, ^°^Mo), with tip distance d = 3.25 fm and 
total excitation energy E* = 2 MeV. 

Fig. 4 Fragments density contour lines for excitation energy E*=2 MeV and 
tip distances d = 2.6, 2.95, 3.10, 3.25 
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Table 1 

Pairs of fragments deformations (/3i,/32) and ratio of excitations energies E2/EI for 

different excitation energies E* and tip distances d. 

E*{MeV) d (fm) (3i (32 E^/El 






2.95 


0.270 


0.370 


- 


2 


2.60 


0.313 


0.492 


0.94 




2.65 


0.325 


0.463 


0.30 






0.263 


0.533 


260. 




2.95 


0.215 


0.486 


0.8 






0.335 


0.351 


0.01 




3.00 


0.209 


0.475 


0.5 






0.332 


0.338 


0.1 




3.10 


0.198 


0.452 


0.2 






0.322 


0.313 


0.5 




3.15 


0.194 


0.438 


0.2 






0.317 


0.302 


1.4 




3.20 


0.309 


0.293 


1.1 


4 


2.60 


0.342 


0.520 


0.70 




3.10 


0.182 


0.530 


0.91 






0.366 


0.324 


0.11 




3.15 


0.175 


0.520 


0.70 






0.362 


0.311 


0.20 




3.20 


0.168 


0.510 


0.53 






0.357 


0.300 


0.30 




3.25 


0.352 


0.286 


0.41 






0.162 


0.499 


0.40 




3.35 


0.341 


0.263 


0.74 






0.151 


0.475 


0.31 




3.40 


0.335 


0.252 


0.98 






0.146 


0.463 


0.14 




3.85 


0.175 


0.265 


0.71 




4.95 


0.335 


-0.212 


0.99 




5.00 


0.337 


-0.226 


0.9 
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Table 1 
(continued) 

^*(MeV) 


d (fm) 


/3i 


(32 


E*2/Et 


6 


2.75 


0.447 


0.409 


0.46 




3.25 


0.158 


0.564 


0.99 






0.388 


0.306 


0.14 




3.30 


0.151 


0.555 


0.81 






0.385 


0.293 


0.21 




3.35 


0.143 


0.546 


0.65 






0.381 


0.279 


0.28 




3.40 


0.134 


0.535 


0.51 






0.377 


0.265 


0.38 




3.45 


0.128 


0.526 


0.43 






0.373 


0.253 


0.48 




3.50 


0.121 


0.516 


0.34 






0.368 


0.239 


0.61 




3.55 


0.114 


0.504 


0.26 






0.364 


0.226 


0.73 




3.60 


0.108 


0.493 


0.20 






0.360 


0.213 


0.87 




5.00 


0.373 


-0.213 


0.49 




5.05 


0.374 


-0.227 


0.46 




5.15 


0.371 


-0.250 


0.52 




5.20 


0.370 


-0.262 


0.58 




5.25 


0.367 


-0.272 


0.65 




5.30 


0.364 


-0.282 


0.75 




5.35 


0.359 


-0.291 


0.90 




5.40 


0.353 


-0.300 


1.1 



12 



148 



Ba 



104 



^ -1180 
^ -1182 



2 

=q -1184 

^ -1186 



-1188 



-1190 



-1192 



-1194 



-1196 



-1198 



-1200 




-0.2 0.2 0.4 



-869 

-870 
-871 
-872 
-873 
-874 
-875 
-876 
-877 
-878 
-879 




P, 



-0.5 -0.25 0.25 0.5 0.75 

P, 




a^Ba ) 



52 



GO. 



0.6 



0.5 



0.4 



0.3 



0.2 



0.1 



-0.1 



-0.2 



-0.3 




TXE = 2 MeV 



20 

15 

10 

5 



-5 

-10 

-15 

-20 
20 

15 

10 



z 



d = 2.6fm 



104., , 148j, 

Mo + Ba 



t z 
A 



-10 
-15 
-20 



d = 3.10 fm 



104., , 148j, 

Mo + Ba 




^2= .313 



P,= .492 




P,= .322 



P,= .313 



t z 

A 



d = 2.95 fm 



104., , 148j, 

Mo + Ba 




t z 

A 



-20 -10 



10 



d = 3.25 fm 



104., , 148j, 

Mo + Ba 



20-20 -10 




^2= .215 



P,= .486 



^2= .296 



^2= .289 
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